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Abstract 

We have obtained the spectral and spatial field correlation functions, Ce(Alu) 
and Ce(Ax), respectively, from the measurement of the microwave field spec- 
trum at a series of points along a line on the output of a random dielectric 
medium. C£j(Ai<j) and Ce(Ax) are shown to be the Fourier transforms, re- 
spectively, of the time of flight distribution, obtained from pulsed measure- 
ments, and of the specific intensity. Unlike (7e(Ao>), the imaginary part of 
Ce(Ax) is shown to vanish as a result of the isotropy of the correlation func- 
tion in the output plane. The complex square of the field correlation function 
gives the short range or C\ contribution to the intensity correlation function 
C. Longer-range contributions to the intensity correlation function are ob- 
tained directly by subtracting C\ from C and are in good agreement with 
theory. 
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The random phasing of multiply-scattered partial waves in disordered media leads to 
intensity fluctuations whose short-range spatial correlation is determined by the square of 
the field correlation function. This field-factorization term is the leading or C\ contribution 
to the cumulant intensity correlation function C = (51 51') ~ C\ — \ (EE 1 *) | 2 . Here 51 is 
the fluctuation of the intensity / from its ensemble average value, 51 = I — (I), E is the 
electromagnetic field, and (. . .) represents the ensemble average. This approximation gives 
the well-known Gaussian statistics of the field JI],^] and the negative exponential distribution 
of the intensity 0-|5[]. The intensity correlation length 5x sets the scale of the speckle spots, 
while the correlation frequency 5v represents the inverse of the spread in transit time between 
the source and the point of observation. Recent advances in speckle statistics have shown 
however that fluctuations in intensity |§-[HJ) total transmission P|TT|-|rT|| and conductance 
|18|^2| are enhanced beyond the predictions of Gaussian statistics, as a result of long and 
infinite-range intensity correlations, the C*2 and C3 contributions to the intensity correlation 
function, respectively. 

In this paper, we present observations of the field correlation function (EE'*) in the 
spectral and spatial domains in measurements of the microwave field transmitted through 
random media. We relate these to their corresponding Fourier transforms, the time of flight 
distribution of photons and the specific intensity, which gives the angular distribution of 
the intensity in the far field. Taking the complex square of the field correlation function 
gives G\ contribution and makes it possible to consider the difference C — C\, which includes 
only terms beyond the field factorization approximation, which are responsible for enhanced 
fluctuations of transmission quantities. Since higher order contributions than C2 are small in 
these samples, C — C\ essentially equals C2. In the past C<i was found from measurements of 



the correlation function of total transmission p4|,p!5| or from fitting the measured intensity 



correlation function using an assumed functional form for G\ and C2 |23| • Here C2 is directly 
separated from the intensity correlation function by subtracting C\ from C. Theoretical 
expressions for C\ and C2 in the frequency domain including absorption and boundary 
conditions are in good agreement with experiments. 

In order to eliminate instrumental distortions of field and intensity spectra, we normalize 
the field by the square root of the average intensity and the intensity by its average value at 
each frequency and position. In addition, measurements are carried out in a frequency range 
in which the change in scattering parameters is small. This allows us to obtain the functional 
form of the relative degree of correlation. We will henceforth denote the normalized intensity 

1 /2 

1 1 (I) by / and the normalized field E / (I) 1 by E. The cumulant correlation function of the 
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normalized intensity will now be denoted by C, which is the sum of three terms p3,|25],[T2],[T3| : 



C — C\ + C 2 4- C3 • 



(1) 



The C\ term is unity for frequency shift Av = and displacement Ax = 0. The spectral field 
correlation function for a plane wave incident upon a slab or upon a quasi-one- dimensional 
sample is given by p6|-[29| : 



C E (Au) = {E{v)E*{v + Av)) 



sinhgofl sinhai/ 
sinh aa 



(2) 



sinh q L' 

where L is the sample length, V = L + 2zb is the effective length corrected for internal 
reflection, z^ = 2£(1 + R)/3(l — R) is the extrapolation length [|30H 3~2"ll, R is the reflection 
coefficient at the boundary averaged over internal incident angles, i is the transport mean 
free path, a = 1/L a is the inverse of the diffusive absorption length L a = ^Dr a , D is the 
diffusion constant, r a is the absorption time, a = 5£/3 is the randomization distance from 
the plane in which the intensity inside the medium extrapolates to zero, go = 7+ — *7-> 
7 2 = K^o 3 + /3 4 ± a 2 ) and (3 = ^2ttAu/D. The C 2 term is of order 1/g at Av = 
and Ax = 0. Far from the localization transition, the expression for the C*2 contribution 



including absorption is given by [j33| ,|3~4 

C 2 (Au) 



4 V" 4 + 7+ sinh (2 7+ L) - 7 _ sin (2 7 _L) - a sinh (2aL) 



(3) 



gL (3 4 cosh (2 7+ L) - cos (2 7 _L) 

where g = Ak 2 £/3irL is the ensemble averaged dimensionless conductance for a tube of 
length L and cross section A. This expression differs from Eq. (14) in by a factor of 
2 to take into account the one-channel-in, one-channel-out measurement described bellow. 
Thus, 



C 2 (Au = 0) 



sinh2aL - 2aL(2 - cosh 2aL) 
AgaL sinh 2 (aL) 



(4) 



In the absence of absorption, C 2 (Ap) = ^. Equations (y) and (f|) do not include inter- 
nal reflection. The complete expression for C 2 including internal reflection is detailed in 
the Appendix A. The C3 contribution is responsible for universal conductance fluctuations 
measured in electronic systems ]r8|-p2T|]. It has been recently observed in optics in time cor- 



relation experiments |35 |. This contribution is of order 1/g 2 and is more strongly suppressed 
by absorption than the C 2 term || and will therefore not be considered here. 

The system studied here is composed of randomly positioned |-inch polystyrene spheres 
at a volume filling fraction of 0.52 contained within a one meter long, 7.6 cm diameter copper 



3 



tube. Wire antennas are used as the emitter and detector at the input and output surfaces 
of the sample. A Hewlett Packard 8722C vector network analyzer is used to measure the 
microwave field, giving its amplitude and phase. Measurements are made in the Ku and K 
bands with frequency steps of 625 kHz and displacements of 1 mm along a 4 cm line running 
symmetrically about the center of the output surface. After the spectrum is taken at a given 
position, the detector is translated by 1 mm. Measurements are made for two positions of 
the input antenna separated by 3 cm. Once spectra are taken at each point on the line, the 
copper tube is rotated briefly to create a new configuration. Measurements are made in 200 
sample realizations and spectral and spatial correlation functions are computed from these 
data. 

The real and imaginary parts of the field correlation function with frequency shift are 
shown in Fig. [j]. Averaging has been performed over sample configuration, position of 
the source and detector and over frequency between 16.8 and 17.8 GHz, where transport 
parameters are known to vary slowly fl3~2| . The real part of the correlation function falls 
quadratically, whereas the imaginary part rises linearly with frequency shift, for small shifts 
From a three parameters fit using Eq. (EI) and the value of the reflection coefficient 



R = 0.13 measured in [32|, we find the diffusion constant D = 3.3x 10 10 cm 2 /s, the absorption 
length L a = 33.3 cm, the penetration depth a = 18.3 cm, which in turn give the extrapolation 
length Zf, = 9.5 cm and the mean free path i = 11.0 cm. 

We have also carried out measurements of pulsed propagation in the same sample. The 
response to a 1 ns pulse shaped by a pulse forming network and mixed to a local oscillator 
at 19 GHz is measured using high bandwidth B&H Electronics amplifiers. The results are 
collected using a digital sampling oscilloscope. The complex square is taken and averaged 
over 4096 configurations and shown in Fig. |2|. Because the incident pulse is much shorter than 
the typical traversal time through the medium, and encompasses a bandwidth much greater 
than the correlation frequency, this measurement gives the time of flight distribution of 
photons propagating through the sample at a carrier frequency of 19 GHz. The line through 
the data is the Fourier transform of the field correlation function which is averaged over 
configuration and frequencies between 18.5 and 19.5 GHz. The agreement between these 
sets of data shows that the field correlation function and the time of flight distribution are 
Fourier transform pairs in agreement with the demonstration in Appendix B. This result is 
analogous to the demonstration that the intensity distribution on the output surface due 
to point excitation (the point spread function) and the intensity correlation function with 
identical shifts in the incident and scattered wave vector [0,0 are Fourier transform pairs 
f3"8fl . In both cases the variables of the incident and outgoing waves are shifted by the same 
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amount. This result holds even in the presence of significant long-range intensity correlation. 
In previous optical studies, |26|| , only the intensity was measured. Its correlation function 
with frequency was shown to be the square of the Fourier transform of the time of flight 
distribution. In that case, long-range correlation was negligible and the measured correlation 
function was essentially equal to C\. 

The real and imaginary parts of the field correlation function with displacement 
Ce(Ax) = (E(x, u)E*(x + Ax, v)) averaged over sample configuration, frequency v between 
18 and 18.5 GHz, and position x, are shown in Fig. |3|. The imaginary part of the field corre- 
lation function with displacement is small, in contrast to the field correlation function with 
frequency (see Fig. p. As shown below, this is a consequence of the isotropy of the random 
field. Such isotropy may be expected at the output of a quasi-one-dimensional sample, far 
from the boundary, or on length scales smaller than the sample thickness at the output of a 
random slab. By the translational invariance of the field correlation function with displace- 
ment, C E (Ax) = (E(x)E*(x + Ax)) = (E(x - Ax)E*(x)) = C%(-Ax). By the isotropy of 
Ce(Ax) with regard to the direction of the displacement, we can exchange — Ax and Ax. 
This gives Ce(Ax) = C* E (Ax). Since Ce(Ax) equals its complex conjugate, its imaginary 
part vanishes. In contrast, the imaginary part of the field correlation function with fre- 
quency shift Ce(Au) does not vanish because the change of Auj to — Au also requires that 
the complex conjugate be taken. 

In transmission, the angular distribution of the scattered intensity is predicted to be |BU 



1(6) = AcosO + cos 2 6 , (5) 

where A = z^/l and the scattering angle 9 is measured from the normal to a reference 
plane £ close to the output surface. Following Ref. [j39j1 , the 2-D Fourier transform of this 



expression is the 2-D spatial correlation function of the field over S. When taken along a 
line, this reduces to 



Ce(Ax) ' 



(1 + 2A) 



^ sm(kpAx) J x (k Ax) 



k()Ax k$Ax 



(6) 



where k is the free space k-vector. This expression gives a good fit over 30 points to the 
data (Fig. |3|) for A = 0.73 and k Q = 3.6 cm -1 , in good agreement with the calculated value 
at 18 GHz of 3.8 cm -1 . This confirms that the scattered intensity angular distribution 1(9) 
and the spatial correlation function Ce(Ax) are Fourier transforms. 

Squaring the field gives the intensity and allows us to compute the intensity correlation 
function and to compare it to the modulus square of the field correlation function. In 
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the frequency domain, we compare the measured intensity correlation function C(Au) to 
C\ + C 2 given in Eqs. (0) and (|A1|), assuming C 3 is negligible and using the value of L a , D, 
a and z b found from the fit to Ce without additional adjustable parameters (Fig. The 
value C(Av = 0) = 1.13 is the variance of the normalized intensity. The difference of this 
from unity is a measure of mesoscopic corrections to correlation. The C — C\ contribution 
normalized by its value at Av = is shown in Fig. |^. The comparison to Eq. ( |A1| ), which 
includes z&, and to Eq. (H) which does not, shows the effect of internal reflection. 

In conclusion, we have measured the field correlation function with frequency and dis- 
placement and shown that they are the Fourier transforms of the time of flight distribution 
and the specific intensity, respectively. Comparison with theoretical expressions gives the 
diffusion coefficient, diffusive absorption length, extrapolation length and mean free path. 
The measurement of both the field and intensity correlation functions makes possible a direct 
separation between short and long-range intensity correlation. 
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APPENDIX A: 

We give the complete expression for the long range contribution Ci including absorption 
and internal reflection, for quasi-one dimensional geometry. We follow the approach of 
in which the correlation function was written in the form: 

m =2 N - + N > + N ' , (A!) 
gBL 

with an overall factor of 2 to account for the one channel-in, one channel-out experimental 
setup. We obtain the following expressions, which differ for B and N a from ||40|| : 
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B 



1 + 2z 2 b ( 1 2 + - 7!) + zt{ll + 7 !) 2 ] [cosh(2 7+ L) - cos(2 7 _L)] 
+ Azl^l + 7!) [cosh(2 7+ L) + cos(2 7 _L)] 

+ [A 1+ z b [l + z 6 2 ( T 2 + 7! )] S inh(2 7+ L)} + {4 7 _z b [l - z 2 ( 7 2 + 7!)] sin(2 7 _L)} , (A2a) 



2 7+ ( 7 ^ - a 2 ) 

x { [(2tJ - a 2 ) + ^ 2 (8 T t + 2 7 ^ 7 ! - lW - ZlW + « 4 ) + 4« + ~ + « 4 ) 

sinh(27 + L) 



+ 



N h 



4z fe 7 + (37i - a 2 ) + 4^ 7 +(37+ + 7^1 - ^iW + « 4 )j sinh 2 ( 7+ L)} 
1 



(A2b) 



2 7 _( 7 2 +a 2 ) 

x { [-(2 7 ! + a 2 ) + z 2 (8 7 l + 27+7- + lW + 3 7 V + a 4 ) - zt^l + 7 !)(2 T i + 2 7 2 a 2 + a 4 )' 
sin(27_L) 



+ 



4z fe 7_(3 7 l + a 2 ) - 4^7_(37 4 „ + 7 2 7 1 + 2 7 V + a 4 ) sin 2 ( 7 _L)} 



(A2c) 



N. 



7+ + 7_ 



" " 2a( 7 2 - a 2 ) (7H. + a 2 ) 
x { [-a 2 - 2 b 2 (67 2 7! + 57 2 a 2 - 57V + a 4 ) + 4a 2 (-2 7 2 7 2 - 7 2 a 2 + 7 2 _a 2 )j sinh(2aL) 



+ 



4z fe a- 1 (-7 2 7l - 7 2 a 2 + 71V - a 4 ) + 44a(-37 2 7i - 2 7 2 a 2 + 27la 2 ) sinh 2 (aL)} . (A2d) 



APPENDIX B: 

We show that the time of flight distribution and the field correlation function with 
frequency shift are Fourier transform pairs. The response to a wave packet f(t) is E(t) = 
J duexp[—i27rut]F(v)E(u), where F(v) is the Fourier transform of f(t) and E(v) is the field 
measured at a given point at frequency v. From this expression, the ensemble average of 
the temporal intensity variation of the transmitted pulse may be written as, 

(E 2 (t)) = J c/z/ 1 rfz/ 2 exp[-z27r(z/ 1 - v 2 )t]{E{v l )E*{v 2 ))F{v l )F*{v 2 ) . (Bl) 

If f(t) is a short pulse with a bandwidth much wider than the correlation frequency Su, F(ui) 
and F(v 2 ) are nearly identical as long as v 2 — V\ is much smaller than the pulse bandwidth. 
With the change of variables v\ = v + Av/2 and v 2 = v — Av/2, F(vi)F*{y 2 ) m \F(u)\ 2 and 
Eq. ( pi| ) becomes 

(E 2 (t)) = [ dudApexp[-i2irAut)(E(p + Av/2)E*{v - Av /2))\F(v)\ 2 . (B2) 
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Within a frequency range in which the variation in the averaged intensity is small, we 

can write (E(u+Au/2)E*(u—Au/2)) = (I(v))Ce(v, Av), where C E (v, Av) is the normalized 
field correlation function around the frequency v. In a frequency range in which the change 
in scattering parameters is small, Ce(v,Av) = Ce(Au) and we find 

(E 2 {t)) = J du(I{u)\F{u)\ 2 ) J dAuexp[-i27iAu t}C E (Au) (B3) 

roo r 

= J dt(I(t)) J dAuexp[-i27rAu t]C E (Au) , 

where the last step follows from Parseval's theorem. Consequently, the time of flight dis- 
tribution T(t) = (E 2 (t))/ J dt{I(t)) and the field correlation function are Fourier transform 
pairs: 

J(t) = / dAv C E {Av)^\-i2ixAv t] . (B4) 
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FIGURES 



FIG. 1. Real (crosses) and imaginary (dots) part of field correlation function with frequency 
shift. Eq. (||) has been used to fit (solid line) the real part. The fit to the imaginary part is not 
shown, but it necessarily follows. 



FIG. 2. Time of flight distribution for carrier frequency 19 GHz. The continuous line is the 
Fourier Transform of the field correlation function (Eq. (B4)) in the range 18.5 to 19.5 GHz 



FIG. 3. Real (circles) and imaginary (dots) part of field correlation function with displacement. 
Eq. (||) has been used to fit (solid line) the real part. The imaginary part is predicted to vanish. 



FIG. 4. Intensity correlation function with frequency shift (circles). The solid line is the theo- 
retical expression for C\ + C2 using the values of L a , D, a and Zb found from the fit of Ce- The 
C\ and C2 contributions are represented by the dotted and dashed lines, respectively. 



FIG. 5. Long range contribution to intensity correlation function C — C\ normalized to its value 
at Au = 0. The solid line corresponds to Eq. ( |A1| ) which includes internal reflection. The dashed 
line is given by Eq. (|3]), which does not include internal reflection. 
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